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Abstract 



Imaging with a spherical mirror in empty space is compared with the case 

when the mirror is fihed with the medium of Maxwell's fish eye. Exact time- 

^^ dependent solutions of Maxwell's equations show that perfect imaging is not 

c/3 achievable with an electrical ideal mirror on its own, but with Maxwell's fish 

. ^H eye in the regime when it implements a curved geometry for full electromag- 

p> I netic waves. 

q 
• ^ 1 Introduction 

Js^ The prospect of perfect imaging with negative refraction [IJ has initiated the en- 

, ^ tire research area of metamaterials that, in turn, inspired the development |2l |3] 

of transformation optics IH El El [7j, the subject of this Special Issue. Yet, ironi- 
►^ cally, negatively-refracting lenses have never perfectly worked in practice, only "poor 

OO man's lenses" , that are substantially thinner than the wavelength, have shown sub- 

Q wavelength imaging |8]. The reason is that negative refraction is only possible in 

O highly dispersive and hence highly dissipative materials |9]; here absorption does 

t^^ not only reduce the intensity but severely limits the resolution of the, theoretically, 

^ perfect lens. Alternatives are hyperlenses [101 HI] that rely on materials with in- 

,-H definite metric. These are anisotropic materials where one of the eigenvalues of the 

kl electric permittivity is negative; these materials thus implement a hyperbolic ge- 

ometry [H] (hence the name hyperlens). Hyperlenses are able to funnel out light 
from near-fields without losing subwavelength detail, but their resolution is given 
by their geometric dimensions, and is not unlimited. 

Another example of perfect imaging has been known, as a theoretical idea, since 
an 1854 paper by Maxwell [12]. This device, called Maxwell's fish eye, because it 
reminded Maxwell of the eyes of fish, uses positive refraction. Maxwell's fish eye 
focuses all light rays emitted from any point in an exact image point; it makes a 
perfect lens for light rays. Luneburg [13] discovered that it maps light rays in physical 
space to rays on a virtual sphere, a curved space; Maxwell's fish eye thus is an early 
example (and inspiration) of non-Euclidean transformation optics [T^. However, it 
is known since Abbe's theory of imaging that the resolution of optical instrument is 



limited by the wave nature of light [To]. Ray optics is not sufficient here, especially in 
curved geometries, unlike in Euclidean transformation optics [3] where geometrical 
optics [15] is exact [HI [7]. Analytic solutions of Maxwell's equations for Maxwell's 
fish eye [161 E] proved theoretically that perfect imaging is possible with positive 
refraction and recent experiments [18j indicate that it also works in practice. 

Yet perfect imaging with positive refraction [THl [IZl Hn] challenges [201 Ell 1221 [23] 
some of the accepted wisdom [T] of subwavelength imaging: it does not appear to 
perform an amplification of evanescent waves [T] and requires a drain at the image in 
the stationary regime [TG]. To illuminate some of the perceived paradoxa of perfect 
imaging, it is instructive to consider the simplest possible case [23j: imaging with 
a spherical mirror. In this paper, the imaging with such a mirror is compared with 
the imaging in Maxwell's fish eye [IT]. Exact time-dependent solutions of Maxwell's 
equations reveal the similarities and characteristic differences between the two cases; 
the mirror cannot perfectly image, but Maxwell's fish eye can. 




Figure 1: Spherical mirror. In this paper we consider the electromagnetic wave emitted 
by a point dipole (double arrow) in the centre of a perfect spherical mirror (half sphere — 
only half of the mirror is shown for being able to see the interior) . We compare this case 
with imaging in Maxwell's fish eye in three dimensions |17j . 



Imagine a point dipole placed at the centre of a spherically symmetric mirror 
that completely surrounds it (Fig. [I]). The dipole emits electromagnetic radiation 
that reaches the mirror, is reflected there, travels back and focuses at the centre 
where the dipole sits — the dipole is imaged to itself. Due to the symmetry of 
the sphere, this is the simplest conceivable imaging problem. Within the regime 
of geometrical optics [12], the focusing with the mirror would be perfect, but the 



highly concentrated field in the centre violates the validity condition of geometrical 
optics [15J. Therefore, exact calculations of the complete Maxwell equations are 
required for elucidating the imaging performance of the mirror. The same applies to 
Maxwell's fish eye [12] where we fill the space enclosed by the mirror with a dielectric 
material that has a certain radially symmetric graded refractive index n{r). In both 
cases, the propagation of electromagnetic waves is governed by a length scale, the 
size of the spherical mirror, and an associated time scale, the travel time of light in 
the mirror. Our results will only depend on these scales. Therefore, to simplify the 
algebra, it is wise to measure space in units of the mirror radius; in these units the 
mirror surrounds the dipole at radius 

r = 1 (mirror). (1) 

It is also advantageous to measure the time t in units of the distance traveled in free 
space; in these units we obtain for the speed of light in vacuum 

c=l. (2) 

Note that in our units the (angular) frequency u is identical to the wavenumber 
in empty space. We consider the propagation of electromagnetic waves in three- 
dimensional space, a case where the theory is mathematically simpler than in two- 
dimensional wave propagation. In 3D, however, impedance matching is required for 
perfect imaging where the electric permittivity e equals the magnetic permeability 
/x, which is difficult to achieve for materials in practice (but trivially the case in 
empty space). For the sake of theoretical simplicity we require 

e = ft = n{r) . (3) 

In empty sphere, the refractive index is unity: 

n = 1 (empty space) (4) 

and in the case when we fill the sphere with Maxwell's fish eye we have 
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n = ^ (Maxwell's fish eye). (5) 

1 + r'^ 

Equations ([T])-(|5]), together with Maxwell's equations of electromagnetism, set the 
scene for the problem we investigate in this paper: the radiation of a dipole placed 
at r = 0. But before we address this problem a few (well-known) ideas on causality 
are needed. 



2 Causality 

Maxwell's electromagnetism is completely time-reversible; a dipole may emit an 
electromagnetic wave, but it may as well absorb the same wave run in reverse 
— by carefully choosing the initial conditions at the boundary an electromagnetic 
wave may focus at the dipole and be completely absorbed. The latter solutions of 




Figure 2: Schematic space-time diagram of retarded waves (left) and advanced waves 
(right). The retarded wave is causal; it is created at some moment in time indicated by 
the planar sheet in the space-time diagram. Prior to that time the wave amplitude is 
strictly zero. On the other hand, the advanced wave converges into a single spatial point 
and vanishes there; it is zero after a certain time. Maxwell's equations have both retarded 
and advanced solutions and superpositions of the two, but only the retarded waves are 
causal. 



Maxwell's equations, where the interaction of the wave with the dipole lies in the 
future, are called advanced, whereas the former solutions, where the dipole causes 
the emission of the wave, are called retarded (Fig. |2|. Clearly, in our case retarded 
solutions are required; they are causal: no radiation is present prior to the event of 
emission at, say, time t = 0. We thus require for all field quantities G: 



G{t) =0 for t < . 
Suppose we represent the quantity G in terms of its Fourier components G: 



Git) 



+ 00 



G{u)e-'^'du 



(6) 



(7) 



The causality condition (|6j) is met if G{u!) is analytic on the upper half complex cu 
plane and decays sufficiently fast, such that we can close the integration contour of 
the Fourier integral ([T]) on the upper half plane for t < when exp(— iwt) decays for 
Imo; > 0; if ^(a;) is analytic the closed-contour integral is zero and G vanishes for 



t < 0, as required (Fig. |3]). Conversely, in the causal case (|6]) we obtain from the 
inverse Fourier transformation 



Giu) 



1 

2^ 



+ 00 



G(t)e^"* dt 



that G{u) obeys the defining Cauchy-Riemann differential equations of analytic 
functions, because the only frequency-dependent term exp(ia;t) in the integral ([s]) 
is analytic, provided the integral ([s]) converges. Convergence is guaranteed, because 
Re(iwt) < for Imw > and t > and so the integrand exponentially decays. 
Furthermore, in the limit Imcu — ?■ oo the term exp(ia;t) tends to zero, so G decays at 
cxD on the upper half plane. We conclude that the Fourier transform of causal fields 
must be analytic and decaying on the upper half plane. We use this requirement to 
deduce the causal solutions for our imaging problem. 




Figure 3: Causality on the complex frequency plane. The Fourier transform of causal 
waves (Fig. [2]) must be analytic and decaying on the upper half complex io plane. The 
singularities (dots) of the Fourier transform lie on the lower half plane or, as in this picture, 
on the real axis. In such a case, the contour (grey line) of the Fourier integral (JTl) can 
be closed on the upper half plane for i < and the integration gives zero, as required by 
causality (pi). The picture shows the first 11 singularities (28) of the Fourier-transformed 
vector potential of dipole radiation in the spherical mirror (Fig. [T]). 



3 Symmetry 

Consider a radially-symmetric impedance-matched medium ([s]) inside the spherical 
mirror at radius ([I| and suppose that a point dipole is placed at the centre that 
emits a flash of light. To describe the electromagnetic wave emitted by the dipole 
we use spherical polar coordinates {r, 9, 0} with z axis aligned to the dipole. We 
assume that the vector potential A of the wave points in the direction of the dipole 



with an amplitude A that only depends on the radius r (our calculations will show 
that this assumption is consistent with Maxwell's equations). The unity vectoijjex 
in z-direction reads in terms of the basis vectors in spherical polars [5], [7] 

sin6' 

Cz = COS Cr Be . (9) 

r 
We thus require 

A' = A{r,t)(cos9,- — ,o]. (10) 



However, the vector potential is rather a one-form Ai [7] (also known as a covariant 
vector) and not the (contravariant) vector A^. We thus lower the index with the 
metric of the spherical polar coordinates [Sj and write down the ansatz for the vector 
potential as 

Ai = A{r, t) (cos 9, -r sin 9, 0) . (11) 



We obtain the magnetic field B from B = 'V x A where for our ansatz (11) we 
calculate the curl in spherical polars [5J and lower the index [5] . We arrive at 

Br = 0, Bg = 0, B^ = B sin^ 9 (12) 

with the magnetic-field amplitude in the plane of emission (for 9 = vr/2) being 

B = -rdrA. (13) 

We use the symbol dr to abbreviate the partial derivative d/dr and we are going 

to use dt for d/dt. The electric field E we obtain from Maxwell's equation in the 

medium ([3| 

R 

Vx dtnE = 6(r)6(t)e, (14) 

n 

with our units ^ and for a point dipole at the centre pointing in z-direction that 

emits a flash of light at t = with normalized intensity. We thus find for the 

non-vanishing electric-field components 

/ 2B 6(r)6(t)\ ^ ^^ f ^r^B r5{r)5{t)\ ^ , , 
\n'^r^ n J \ n n n J 

while the azimuthal component dtE^ is zero. Furthermore, we use Faraday's law 
^ X E = —dfB expressed in spherical polars [S] with indices lowered [S] and our 



result (12) for the magnetic field, to obtain the wave equation 



n n n'^r'^ n 



All other components of V x dtE vanish for the electric field (15) in spherical 



polars, which reduces the problem to the wave equation (16) and justifies our ansatz 



-'^Note that we employ coordinate bases [Sl[7] and not the frequently used orthonormal systems 
of basis vectors that, apart from Cartesian systems, are non-coordinate bases. Coordinate bases 
are much simpler to work with, given a little knowledge of differential geometry [51 [J • 



(11). Finally, we require that the mirror acts as an ideal electrical mirror where the 



electric-field components in the mirror are put to zero, in our case (jTl) 



Ea 



r=l 



which implies from expression (15) that 

B 



dr 



n 



0. 



(17) 



(1^ 



r=l 



We see that the magnetic field must not change at the mirror. In the following we 



seek causal solutions of the wave equation (16) with boundary condition (18) for the 



magnetic field B expressed in terms (13) of the vector potential. First we consider 
the dipole radiation in empty space surrounded by the spherical mirror j23] and then 
we compare this case with Maxwell's fish eye [17]. 



4 Mirror 



In empty space ^ the left-hand side of the wave equation ( 16 ) for the magnetic 
field (13) reduces to 



dlB - ^ - d^B = -r dr (dlA +-drA- dfA 



(19) 



such that the wave equation ( 16 ) is satisfied if 



d^^A + -drA- dlA = -6{r)6(t) 



(20) 



We thus require that the Fourier components A of the vector potential, with the 
convention (ItI) for the Fourier transform, obey the inhomogeneous wave equation 



d'^A + -drA + 00^ A 



2tx 



(21) 



For r > the general solution of the — then homogeneous — wave equation (21 ) is 

A 



A 



A, 



(22) 



Note again that in our units ^ the frequency co is equal to the wave number that 
normally appears in the spherical waves (22). The first term of expression (22) 



describes an outgoing wave and the second term an ingoing wave. The spherical 
mirror turns outgoing into ingoing waves and so we expect that the coefficients A± 



are not independent. Indeed, we obtain from condition (18) for the magnetic field 



(13) of the wave (22) 



U"^ — lOJ — 1 

7 



-^+e2-+2'^ 



(23) 



where 6 denotes the phase 



S = arctan 



OJ 



w" 



(24) 



At the mirror ([I]) the outgoing component ^+exp(ic<;r) thus produces the ingoing 
component ^_ exp(— iwr) and vice versa, apart from the extra phase (24). This 



phase vanishes for high frequencies a; — )■ oo where the wavelength approaches zero 
and imaging becomes perfect, regardless whether it is subwavelength-limited or not. 



For finite u the phase (24) is responsible for limiting the resolution of the mirror, as 
we see next. 

So far we considered the general solution of the homogeneous radial wave equa- 



tion with the boundary condition (18) at the mirror. Now we turn to the solution 



that satisfies the inhomogeneous equation (21) and therefore describes the field gen- 



erated by the point dipole. For this we write down the following combination of 



outgoing and ingoing waves ( 22 ) 



A 



sin(u; — ur + 6) 
Svr^r sin(co' + 6) 



(25) 



One verifies that this A obeys the condition (23) of spherical waves reflected by the 
mirror. We also see that 

for r ^ , (26) 



A 



^TT^r 



where [Anr) ^ is the Green function of the Poisson equation. 



r 



47rr 



Anr 



-6{r) . 



(27) 



In the vicinity of the origin we can ignore features of A on the scale of the wavelength 



that depend on the frequency oj and so the inhomogeneous wave equation (21) re- 



duces to the Poisson equation here. Therefore the wave (25) with the asymptotics 



(26) satisfies the wave equation (21), and we already know that it also obeys the 



boundary condition (18). The Fourier components (25) thus constitute the electro- 



magnetic field of a light flash emitted by a point dipole, provided they are causal, 
i.e. analytic and decaying on the upper half u plane. 

Causality is the final point we need to consider for the field (25). Representing 



the sine functions in expression (25) in terms of exponentials we see that A — )■ 
53t<j — 7- oo. Furthermore, A is analytic, apart from poles Um on the real axis (Fig 
where 

+ Sm) = 0, Sm = 5{uJrn) ■ (2 



or 



sm Wr; 



The Fourier transform A also has a quadratic singularity at w = 0. In order to 
obtain a causal solution we must move the poles slightly below the real axis, by 
providing u with an infinitesimal, positive imaginary part. Alternatively, we move 
the integration contour up from the real axis by an infinitesimal distance (Fig. pi). 
In both cases, we can close the contour of the Fourier integral ^ on the upper half 



plane for t < and get zero, as required for causal waves (|6j). For t > we close the 
contour on the lower half plane and obtain from Cauchy's theorem 



A = —^ ^ rim sin(wmr) sin(a;mt) 



m=l 



with the coefficients 



UJrr, — (^m + 1 



U. 



4 - 2a;2 



(29) 



(30) 



For the light flash emitted by the dipole in the centre of the mirror, we obtained the 



exact solution (29) with the pole frequencies (28) and the coefficients (30). Let us 



investigate the extent to which this wave images the dipole onto itself. 

For perfect imaging we require that an outgoing wave is perfectly converted into 
an ingoing wave. The wave is focused at the centre and bounces back and forth 
with a period of twice the time it takes to travel to the mirror, 2 in our units. As a 
flash of light is described by a radial delta function, the flash in a perfectly imaging 
device should be given in terms of the periodic delta function A(t) of period 2, i.e. 
the kernel of the discrete Fourier transform 



A(t) 



+ 00 

E 

m,=— oo 



-. +00 ^ oo 

6{t -2m) = - Y, e'""* = n + Yl cos(m7rt) . 



(31) 



171=1 



At the mirror ([I]), the wave changes sign and propagates backwards. We thus require 
for the radial wave Aq in the case of perfect imaging of the dipole onto itself 



^0 



A(t 



A(t 



2) A(t-r) A(t + r) 



(32) 



Arc'^r An^r An^r Air^r 

Writing Aq in terms of the periodic delta function A guarantees that Aq is periodic 



with period 2. We obtain from expression (31) 



An 



2n^r 



y sin(r?7,7rr) sin(m7rt) 



(33) 



m=l 



Formula (33) strongly resembles our result (29) for the spherical mirror, but not 



perfectly. There, instead of the eigenfrequencies rmr the pole frequencies Um appear 



poles (28) 



and the trigonometric series (29) contains the coefficients rjm- The poles (28) depend 



on the phase (24) that, for large u, approaches u ^. In this limit we obtain for the 



UJr, 



± 



mn 



+ 



m?T{'^ 



1 



(integer m) 



(34) 



2 V 4 
that approaches rnvr for m — ?■ oo. In this limit also r^m — ^ 1, as we see from expression 



(30). Consequently, whenever large-m components dominate the series (33) of the 



perfect wave, they also dominate the wave (29) in the spherical mirror. We see from 



definition (32) that this is the case at the infinitely peaked front of the light flash 



bouncing back and forth in the mirror. Therefore, the perfectly imaging wave makes 
the dominant contribution to the electromagnetic wave in the spherical mirror. 



Figure 4 shows the difference A — Aq between the vector potential (29) of the light 



reflected in the spherical mirror and perfectly reflected light flashes (32). Relativistic 



causality implies that the outgoing wave A must remain a perfect pulse until it 
hits the mirror. Therefore the deviation from the ideal vector potential Aq cannot 
correspond to any physical field. Figure ^ indicates that, during this stage. A — Aq 
is constant in space but grows in time. This growth must be linear in order for 



A to obey the homogeneous wave equation (20). As the magnetic field (13) and 



the electric field ( 15 ) depend on the spatial derivative of A in the gauge relevant to 



our ansatz (11), the electromagnetic fields do not deviate from free dipole radiation 



until the wave hits the mirror. After the first reflection the wave is no longer a 
perfect pulse: the vector potential carries an additional discontinuity at the position 
of the pulse and develops a field in its wake. Note that the mirror is not assumed 
to be dispersive or otherwise imperfect; these features quantify the imperfections of 
imaging with an ideal spherical mirror. 
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Figure 4: The reflection of a light pulse in an ideal electrical mirror is not perfect. The 
figure shows AirlA — Aq) where A is the amplitude (29) of the vector potential of dipole 
radiation in the spherical mirror and Aq is the amplitude (32) of a perfectly reflected 



inflnitely peaked pulse. The different shades of grey indicate the different evolution times 
t = 0.25,0.5,0.75,1.0,1.25,1.5,1.75 in our units, the lighter the tone the earlier is the 
time. One sees that before the pulse has reached the mirror at r = 1 the amplitude of the 
radial vector potential is uniform in space but grows linearly in time. After the interaction 
with the mirror the amplitude develops an additional discontinuity at the position of the 
pulse and a fleld in its wake (see also Fig. [5|. 
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5 Perfection 

Now imagine we fill the space enclosed by the mirror with the impedance-matched 
([3]) medium ([s]) of Maxwell's fish eye 112]. This device implements the geometry of 
a curved space. In particular, it creates the illusion that light propagates on the 
3-dimensional surface of the 4-dimensional hypersphere [13] . It turns out |17] that 
the geometry of waves in Maxwell's fish eye appears in the clearest possible form if 
we write the vector potential A as 

A = 2 (drD) n dr or, equivalently, drA = 2ndrD (35) 



or, from definition (13), 



B = -2nrdrD . (36) 

We show next that the function D is the scalar Green function on the hypersphere 



[T7] . For this we inspect the left-hand side of the wave equation (16) where 



dr-dr-- ^ - d^B = ~2nrdr \ -^ dr-nr^drD - D - d^D ] . (37) 



Therefore the wave equation (16) for the refractive-index profile (pi) with n(0) = 2 
is satisfied if D obeys 

^ drnr'drD-D-d^D = -^-^^'^, (38) 






or, in Cartesian coordinates. 



^V.nVD-D-d?D = J-^^. (39) 



w n-' 



This is the conformally-coupled radial scalar wave equation on the virtual space 
implemented by Maxwell's fish eye, the surface of the hypersphere [T7]. If D satisfies 



the boundary condition (18) at the mirror in addition to the wave equation (38) -D 
is the required scalar Green function. 

Let us construct the Green function D. For this we Fourier-transform the wave 



equation (39) 



\v-nVD + {u^-l)D = -f^^, (40) 



n^ 2Txn 

and obtain for the Fourier-transformed radial wave equation: 



n-^r 



drnr'drD + {uj''-l)D = -^^. (41) 



3^2 2T:n^ 



We consider the limit r — )■ where the far-field terms (o;^ — 1)D do not influence 
the near field of D anymore, such that we can replace the left-hand side of the 



Fourier-transformed radial wave equation (41) by 
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if D obeys the asymptotics 

One easily verifies the following expression for the two fundamental solutions of the 



Fourier-transformed wave equation (41) with asymptotics (43): 

1 



D 



± 



r -\ — ) exp(±2iu;arctanr) . 



(47r)2 
Let us write down the expression 

1\ sin(2warccotr) — sin(2co'arctanr) 



(44) 



D 



r + 



(45) 



(47r)2 sin('7ra;) 
We see from Euler's formula exp(ix) = cos x + i sin x and 2 arccot r = vr — 2 arctan r 



that expression (45) is a linear combination of the two fundamental solutions (44) 



and hence a solution of the homogeneous wave equation (41 ) for r 7^ 0. Additionally 



since 2 arccot r — )■ tt for r — )■ cxd expression (45) obeys the required asymptotics (43) 



for r — ^ 0. Furthermore, one verifies that D conforms to the boundary condition 



( 18 ) at the mirror, 



dr-rdr-D 



0. 



r=l 



(46) 



Formula (45) thus qualifies as the Fourier transform of a Green function. Finally, 



D decays for '^u — > ±00 and is analytic, apart from single poles at 

(non-zero integer m) 



UJrt 



m 



(47) 



and a double pole at w = 0. Similar to the case of the empty spherical mirror 
(Sec. |4]), we move the poises slightly below the real axis or move the contour of the 
Fourier integral ([T]) above the real axis (Fig. |3| in order to obtain the causal Green 
function D. 

The causal Green function describes the time evolution of a flash of light emitted 
at the point source in the centre of the flsh-eye mirror. We obtain from Cauchy's 



theorem for the Fourier transform of expression ( 45 ) 
0(t) 



D 



A-K^' 



r + 



y [sin(2m arctan r) — sin(2?7i arccot r)] sm{mt) (48) 



m=l 



and write this result in terms of the periodic delta functions (31) similar to formula 



( 32 ) for the ideal imaging wave 

e(t) / 1 



D 



87r2 



A 
-A 



t — 2 arctan r 



7r 



A 



t + 2 arctan r 



TT 



t — 2 arccot r\ /t + 2 arccot r 

+ A 



vr 



TT 



(49) 
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We write A in terms of the defining series (31) of delta functions, extract a factor 
of 77 and combine in the resuhing series the numbers 2m and 2m + 1 from the A 
series (ISll) in one summation index m,'. We obtain the resuh 



D = —^ ('" + ") X] [^*^^ ~ ^ arctan r - m'vr) -6it + 2 arctan r - m'n)] . (50) 



Formula (50) shows that the light flash emitted by the point dipole propagates 
through Maxwell's fish eye as a perfect flash. It is reflected at the mirror where 
the amphtude changes sign, it focuses at the centre and bounces back and forth 
without disintegrating in time: in contrast to the spherical mirror. Maxwell's fish 
eye perfectly images in time (Fig. Isl). 

6 Imaging 

We have seen that Maxwell's fish eye faithfully transmits time-dependent signals, 
it perfectly images in time, but how does it image in space? The spherical mirror 
focuses only one point, the centre, to itself; an elliptical mirror images one focal point 
to the other focal point. Maxwell's fish eye is a so-called perfect optical instrument 
[T5] where all light rays from any source point Vq focus at the corresponding image 
point Tq. The fish eye is a perfect instrument for light rays, but the optical resolution 
of imaging is normally limited by the wave nature of light. Remarkably, one can 
derive exact solutions of Maxwell's equations for electromagnetic waves in Maxwell's 
fish eye [IT] (that needs to be impedance-matched in 3D) but the technicalities 
[7] involved in justifying and explaining them go beyond the scope of this article. 
The principal technical problem is the tensor nature of the electromagnetic field. 
Fortunately, we can liberate ourselves from the technicalities of tensor analysis and 
understand the essence of perfect imaging with positive refraction by discussing the 
scalar Green function D; the electromagnetic field can be derived from D by certain 



differential forms [ZlIIZl that reduce to expressions (36) and (15) for the special case 
of the source being at the centre of Maxwell's fish eye. For simplicity, let us also 
dispense with the mirror and consider the infinitely extended profile ([s]) for r from 
to oo. In this case we write down the Fourier-transform of the scalar Green function 



similar to expression (45), except that we omit the contribution due to the reflection 
at the mirror, 

n-fr'\ M sin(2a;arccotrO 

^" V rV (47r)2sin(7ra;) ' ^^^' 

and replace the radius r by the Mobius-transformed radius [3 [17] 

r' = '^~^°' = . (52) 

a/1 + 2r • ro + |rP|roP 



One can show that formulae ( 51 ) and ( 52 ) describe solutions of the Fourier-transformed 



wave equation (40) with the correct asymptotics at the source point. 



^ ~ 7T^ 1 • (53) 

(47r)^|r — ro| 
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In the ray optics of Maxwell's fish eye [7] , source and image correspond to the and 



the oo of r', the Mobius-transformed radius (52): 



r' = O r = Vq (source), 

r' = oo 4-7- r = r'n = —- — — (image). (54) 

However, the wave D is not singular at the image, because we obtain from formula 



(51) 



~ ujsmc(2uj/r) . . smx , ^, 

U ~ — with smcx = . (55) 

Stt^ sm(7r(jj) x 

We can interpret D in two ways, as the Fourier transform of the time-dependent 
Green function D that describes a flash of light or, alternatively, as the amplitude of 
the continuous wave generated by a stationary point source that emits and absorbs 



radiation in a stationary state. Our result (55) shows that the electromagnetic wave 
of a stationary source develops a diffraction-limited image, in agreement with exper- 
iment p[8] . However, let us adopt the alternative interpretation of D as the Fourier 
transform of a light flash D. We obtain in complete analogy to the calculations in 
the previous section: 

L) = — ^(r'+-) ^ [(5(t-2arctanr'-2m7r)-(5(t + 2arctanr'-2m7r)] . (56) 

^ ^ m — — rv^ 



m=~cx) 



The flash perfectly focuses at the image point, is reflected there and bounces back 
to the source etc. In Maxwell's fish eye the image point thus acts like a perfect 
mirror. In the stationary regime we can imagine the wave as a continuous stream 
of light flashes. Each of the elementary flashes is reflected at the image point where 
its amplitude changes sign. In the stationary regime we average over the stream 
of flashes and so the sign change upon reflection causes the image to get blurred. 
Maxwell's flsh eye has the potential of perfect imaging, but this potential is not 
realized yet. 

What is missing is a crucial ingredient of imaging: a detector. A detector ex- 
tracts the fleld at the image point. An ideal point detector acts as a completely 
passive outlet with point-like resolution (inflnitely small cross section). Such a de- 
tector could be part of a CCD array or represent the photosensitive molecule of a 
photographic material that happens to be at the image point. Suppose the detector 
is positioned at the point Vq and extracts the radiation incident at that point. In 



this case the series ( 56 ) of flashes bouncing back and forth reduces to one flash that 



is emitted at the source and disappears at the image at time t = vr in our units: 

D = ^(r'+-] 6(t - 2 arctan r') e(7r - t) . (57) 

on \ r' J 

For this expression we obtain the Fourier transform 



= ^2 ( '"' "I — ~,] exp(2iu; arctan r') . (5^ 
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Note that formula (58) describes a running wave with complex wave function prop- 



agating from the source to the image where the wave disappears, in contrast to 



the standing wave (51) with real wave function that is reflected at the image. The 
spatial singularity of D for r' — )■ oo corresponds to a supplementary source at the 
image point Tq, a drain. As 2arctanr' — )■ vr for r' — )■ oo the wave carries a phase 
delay of t^uj at the image ^\. The detector, acting as a drain, creates a perfect 
spatial image of the source point. So far, we assumed that the detector is positioned 
at the correct place. Now imagine the detector is elsewhere. As the detector has 
point-like resolution it can only detect a field that is infinitely concentrated at that 
point. Therefore, an ideal point detector placed at the wrong position will not de- 
tect anything; the perfectly focused field will only appear at the perfectly positioned 
point detectoin in agreement with experiments [IB] . Of course, the detector may 
have a non-perfect efficiency such that it does not extract the entire field. In this 
case, part of the field is reflected, but the transmitted part is perfectly focused; the 
detected image is infinitely sharp. As any continuous distribution of sources can be 
thought of as a distribution of source points, any source is perfectly imaged, as long 
as it is detected by a continuous array of point detectors. Only the detected field is 
imaged with point-like precision, but detection is the very point of imaging. 

However, one may still argue that the drain at the image is an artifact that per- 
turbs the field, creating an infinitely focused field entirely on its own, the illusion of 
perfect imaging but not a true image [201 1221 123] . Let us consider a counter example 
where the introduction of a drain turns out not to improve imaging in time: the 
spherical mirror. Returning to the starting point of this article, we consider empty 
space enclosed by the spherical mirror. We form the following linear combination of 



the two fundamental solutions (22) with the reflection condition (23): 



A=Ar (e'"' - e2"^-*-2'^-"^0 , e^'^ = "^l^'"^ ^ . (59) 

Svr^r uj^ — iuj — 1 

This expression creates the singularity of the source at the centre and a drain with 
phase delay 2uj + 25 in our units for a wave reflected at the mirror with phase shift 



25 back to the source in time t = 2. However, expression (59) also develops poles in 



the complex frequency plane at the zeros oi u'^ — iu — 1, i.e. at 



cj = -^ — , (60) 



on the upper half plane! Therefore, the solution (59) is not causal, the drain is 
not causally connected to the source, but rather acts as an independent source of 
radiation. However, as we have seen. Maxwell's fish eye makes all the difference 
here; in this case the drain is consistent with causality. 



^In practice [THj the resolution is limited by the finite electromagnetic cross section of detectors. 
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7 Credo 

An ideal spherical mirror cannot perfectly image on its own, the mirror inevitably 
distorts light pulses (left of Fig. Isl). However, if the space enclosed by the mirror is 
filled with the medium of the impedance-matched Maxwell fish eye, pulses are no 
longer distorted and thus the imaging is perfect in time (right of Fig. Isl). We have 
seen that the introduction of a drain at the image is allowed by causality and leads 
to perfect imaging in space in the stationary regime (similar to the perfect focusing 
of ultrasound waves by time reversal [23 )• 0^ the other hand, in imperfect imaging 
the drain would be in conflict with causality. One of the puzzles of perfect imaging 
with positive refraction seems to be resolved. 

Yet one may still wonder how Maxwell's fish eye is able to restore evanescent 
waves in a similar way negative refraction does [1]. But note that the amplification 
of evanescent waves is not the only physical picture explaining the performance of 
the negatively-refracting perfect lens; there is a simple geometrical argument [25] : 
the lens performs a folding transformation of space [B]. The negative-index lens 
appears to fold the space perceived by electromagnetic waves such that three re- 
gions are identical, one region in front of the device, the lens itself, and the region 
behind it (Fig. [6]). As these regions are identical for the field, the field strength of an 
electromagnetic waves must be exactly the same at all triples of connected points: 
negative refraction makes a perfect lens by spatial transformation [2S]- In this geo- 
metrical picture of imaging we do not need to discriminate between evanescent and 
propagating waves, because the electromagnetic field is transformed in its entirety. 

Geometry also explains the perfect imaging with Maxwell's fish eye that has 
a positive refractive-index profile. Like in the case of the perfect lens, this device 
changes the geometry of space for the electromagnetic field. It creates the effect that 
electromagnetic waves propagate in the 3-dimensional surface of the 4-dimensional 
hypersphere. Imagine an ordinary 3-dimensional sphere with 2-dimensional curved 
surface (the hypersphere is not much different). Suppose light is confined to the 
surface. There it is completely natural that light waves emitted at any point of the 
sphere converge on its antipodal point (left of Fig. [T]). In Maxwell's fish eye, the 
antipodal point on the virtual sphere corresponds to the image point in physical 
space. In the case discussed in this paper, we cut the fish eye by a spherical mirror 
that, in virtual space, appears as a mirror around the equator of the sphere (right of 
Fig. [T]). The introduction of the mirror is not essential for perfect imaging in theory, 
but it is vital in practice, for limiting the index range and avoiding superluminal 
propagation [TB]. The mirror creates the illusion that light explores the entire sphere, 
whereas in reality it is confined to the semi-sphere where the index range in physical 
space is smallest. Such conjuring tricks with mirrors are only possible in curved 
space. Moreover, the imaging of Maxwell's fish-eye, with or without mirror, is not 
restricted to a single point but happens for all points. This is possible for devices 
that implement perfect, curved geometries beyond the approximation of geometrical 
optics. In such cases, the virtual geometry of light created by the material is valid 
for full electromagnetic waves. Perfect imaging is not done with mirrors [26], but 
with geometry. 
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Figure 5: Mirror versus Maxwell's fish eye. Space-time diagrams of the vector poten- 
tial of light flashes evolving in time (left) for imaging in a spherical mirror and (right) in 
Maxwell's fish eye supplemented by a mirror. The mirror is placed at the vertical bound- 
aries of each diagram where the grey vertical line indicates the world line of the centre. 
The black curves correspond to the infinite peaks of the flashes, the contours describe 
the finite values of the vector potential. As the right picture contains no finite contours, 
just infinite peaks, the flash is not disintegrating and bounces back and forth between the 
mirror in Maxwell's fish eye: imaging is perfect there. In contrast, the left picture shows 
the imperfections arising from reflection in an ideal electrical mirror: after the flash has 
bounced off the mirror the vector potential carries an additional discontinuity and a field 
in its wake that develops a singularity at the refocusing time i = 2 in our units (the con- 
tours are cut off around the singularity) . After that time the entire interior of the mirror 
is filled with a spatially non-uniform vector potential that corresponds to a non-vanishing 
electromagnetic field distorting the flash: perfect imaging is not done with mirrors. 
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Figure 6: Negative refraction and folding of space. The figure illustrates the geometrical 
representation [25] of perfect imaging by negative refraction: the device appears to fold 
space such that the entire electromagnetic field in the folded regions is identical — the 
device transfers the field from one region into another: imaging is perfect. Note that this 
geometrical picture does not discriminate between propagating and evanescent waves. 





Figure 7: Geometry behind Maxwell's fish eye. Luneburg |13j discovered that Maxwell's 
fish eye implements in 2D the geometry of the two-dimensional surface of a 3D sphere 
(or in 3D the three-dimensional surface of the 4D hypersphere). The figure illustrates 
the 2D case |16j . The left picture shows a wave on the sphere emitted from a point 
source and forming a perfect image at its antipodal point. In the right picture a circular 
mirror is placed at the equator of the sphere such that the image remains on the southern 
hemisphere [TE] that corresponds to a finite region in physical space. In 3D a spherical 
mirror should be used. Source and image merge when the source is placed at the south 
pole, which in physical space corresponds to emission from the centre of Maxwell's fish 
eye, the case considered in this paper. 
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